Study of energy transfer in helium atom scattering from surfaces 
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Recently developed quantum mechanical theory of inelas- 
tic He atom scattering (HAS) from solid surfaces is employed 
to analyze the energy transfer between projectile particles 
(thermal energy He-atoms) and vibrational degrees of free- 
dom (phonons) characteristic of a variety of experimentally 
studied surfaces. We have first calculated the angular re- 
solved energy transfer which can be directly compared with 
the values deducible from the HAS time-of-flight spectra and 
a good agreement with experimental data has been found. 
This enabled us to calculate the total or angular integrated en- 
ergy transfer, which is of paramount importance in the stud- 
ies of gas-surface scattering, but is neither accessible in HAS 
(which yields only the angular resolved quantities) , nor in the 
wind tunnel measurements for surfaces whose atomic com- 
position and cleanliness must be maintained during the ex- 
periment. Here we present the results for prototype collision 
systems of this kind, viz. He^Cu(OOl), He ^Xe/Cu(lll) 
and He— >Xe(lll) which are representative of the very differ- 
ent types of surface vibrational dynamics and thereby provide 
an insight into some common properties of energy transfer in 
gas-surface scattering. 

PACS numbers: 68.35. Ja, 34.50.Dy, 63.22.-|-m, 47.45Nd 



I. INTRODUCTION 

For a long time the energy transfer in gas-surface col- 
lisions has been in the focus of interest of surface scien- 
tists. Special attention has been paid to the energy (heat) 
transfer in the free molecular flow regime because of its 
importance for understanding the flow interaction with 
space vehicles. The energy transfer in this regime of gas- 
surface collisions is dominantly due to the gas interaction 
with vibrational degrees of freedom or phonons of the sur- 
face involved. The macroscopic features of heat transfer 
characteristic of the various technologically interesting 
surfaces are routinely investigated in wind tunnel exper- 
iments but the microscopic properties of the correspond- 
ing interactions can be assessed only by molecular beam 
scattering techniques. 

The major amount of information on the microscopic 
vibrational properties of a number of surfaces of both 
fundamental and practical interest has been acquired by 
utilizing the technique of mono-energetic He atom scat- 
tering (HAS) combined with the energy resolution of the 
scattered beam by the time-of-flight (TOF) method (for 
a review see Ref. [ |l|]). At present, this technique pro- 
vides only the angular resolved quantities (e.g. energy 



transfer in angular resolved scattering events), and in 
order to obtain experimental information on the total 
energy transfer in gas-surface collisions the HAS studies 
must be complemented by wind tunnel measurementsn. 
However, the studies of total energy transfer at surfaces 
whose structure, composition and cleanliness can be con- 
trolled only under the UHV conditions, are not possible 
in the wind tunnel environment. Hence, different meth- 
ods for assessing the total energy transfer involving such 
surfaces must be employed in order to obtain the relevant 
information. 

In this work we demonstrate that investigations of 
vibrational properties of surfaces by the application of 
HAS-TOF technique also provide an excellent database 
for the assessment of total energy transfer in gas-surface 
scattering. Whereas the experimental data acquired in 
the highly quantum single phonon collision regime of 
HAS are primarily used to detect and identify the var- 
ious surface localized vibrational modes, the available 
information from the single and multiphonon scattering 
regime combined with the quantum theory of HAS can be 
conveniently used to reveal the magnitudes of angular re- 
solved and total energy transfer in gas-surface collisions. 
This is particularly important for a relatively large num- 
ber of surfaces whose desired microscopic characteristics 
cannot be maintained in wind tunnel experiments which 
at present are the only method that can render the values 
of total energy transfer. 

In the following sections we first demonstrate how the 
quantum, microscopic theory of energy transfer in gas- 
surface collisions can be verified by comparing with the 
experimental data obtained from HAS-TOF measure- 
ments. Then we show how such a theory can be extended 
to calculations of the total energy exchange between the 
projectile particles and the target as a function of the sur- 
face temperature, projectile incident energy and angle, 
etc. The procedure is exemplified on selected prototype 
systems recently studied experimentally: clean Cu(OOl) 
surface, monolayer of Xe atoms adsorbed on Cu(lll) sur- 
face and the (111) surface of Xe multilayer crystal con- 
densed on Ft (111). The properties of energy transfer 
involving these surfaces with manifestly different types 
of vibrational dynamics should be equally relevant to a 
variety of other systems investigated by HAS, as well as 
to the surfaces studied in wind tunnel experiments. 
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II. QUANTUM DESCRIPTION OF 
ATOM-SURFACE SCATTERING AND ENERGY 
TRANSFER 

The knowledge of the quantum scattering spectrum 
Mci.Ts(£) which gives the probability density that an 
atom with initial momentum Tiki exchanges energy e with 
a surface at the temperature Ts can be used to evaluate 
the total energy transfer fi, wbicjii enters heat transfer 
and accommodation coefficientso according to the ex- 
pression 



(1) 



The quantity Nk^^T^{e), and thereby /i(ki,rs), is not 
directly accessible in typical HAS-TOF measurements 
from which most of the data is available at present. 
These experiments yield the energy and angular resolved 
quantities usually only for fixed total scattering angle 
SsD = Si + 9 j where 9i and 9j denote the incident and fi- 
nal angle of scattering, respectively. The TOF spectrum 
is directly proportional to the energy and lateral momen- 
tum resolved scattering distribution iVki.T^ (e, AK) (see 
Ref. [ H ) in which e and AK denote the amount of energy 
and lateral momentum exchanged between the particle 
and the surface. The values of e and AK are connected 
through the conservation of total lateral momentum and 
energy in the course of collision, leaving e and Of as the 
only independent observation variables. Using this we 
can obtain the angular resolved energy transfer from: 



/eiVki,T,(£,AK(£))d£ 

/7Vki,T.(£,AK(£))d£ ' 



(2) 



and the magnitude of /ir(ki, T^, 0/) computed from the 
experimental TOF spectra can be compared with the re- 
sults based on theoretical A'kj.Ts (£, AK). The calcula- 
tion of reliable multiphonon A'ki,Ta (e, AK), from which 
-^ki,Ts(£) is obtained by integration over AK, is dif- 
ficult in the regimes in which multiphonon and quan- 
tum dynamics effects are important. However, recent 
progress in experiment and theory of HAS from surfaces 
makes it possible to measure and accurately calculate 
iVkj,Ts (e, AK) for some prototype collision systems and 
thus obtain a reliable estimate of the total energy transfer 
in these experiments. 

In the regime of HAS in which the uncorrelated phonon 
processes are dominant the angular res(a||Yed scattering 
spectrum has the following unitary forr 



. AT^N f drd 
iVk,T.(e,AK)= / — 



drd^R 



-oo (27r?i)3 

X eH^r-h{AK.)Ti.) exp[2T4^(T, R) - 2W]. (3) 

Here 2W{t, R) is the scattering function to be defined bn^ 
low, 2W = 2VF(0, 0) gives the Debye- Waller exponent^ 
and T and R are auxiliary integration variables used to 



project the states with e and AK out of the integral on 
the RHS of (|3|). In the range of validity of (|) the en- 
ergy and lateral momentum are conserved in each phonon 
exchange process and the values of e and AK are con- 
strained by the conservation of the total lateral momen- 
tum and energy in the collision. Using expression (^) it 
is possible to write Eq. p) as: 



d 

i—2W(T 

OT 



0,R = 0), 



(4) 



which can be readily calculated once 2W(t, R) is estab- 
lished. In the following the projectile-phonon coupling is 
assumed linear in the adsorbate displacements since this 
gives the dominant multiphonon contribution observed 
in HAsIa, yieldingQ: 



2M^(r,R)= ^ [\Vf:fJ{+)\'[n{cu^,) + l]e 



-i(wQjT — QR) 



V, 



(-) p n(cjQj) 



,i(wQjr-QR) 



(5) 



Here Q and j denote the lateral wave-vector and branch 
index of a normal phonon mode of frequency ujQj , respec- 
tively, k ~ (K, fc^) where kz is the quantum number of 
distorted waves describing the projectile motion normal 
to the surface and n{u!Qj) is the Bose-Einstein distribu- 
tion. The symbols | V^';*^:-'(±) p denote one phonon 
emission (-I-) and absorption (— ) scattering probabilities 
expressed in terms of the corresponding on-the-energy- 
shell transition matrix elements V^'^'"^.'"' (±_) normalized 
to unit particle current in the z-directionQ. Substitu- 
tion of (||) into (|^) gives the desired expression for the 
mean energy transfer in which the so-called recoil term 



determines the temperature dependence. Quite gener- 
ally, the magnitude of the total energy transfer calculated 
from expression is much more sensitive to the char- 
acter and variations of the polarization vectors of surface 
modes than to the changes in the functional behavior of 
their dispersion. A more detailed description of the scat- 
tering formalism used in the present study has been given 
in Ref. [§. 

It should be pointed out that the above quoted expres- 
sions for calculating the scattering spectra and energy 
transfer are valid irrespective of the collision regime (sin- 
gle versus multiphonon) they are applied to. Hence, the 
reliability of the scattering potentials and the description 
of the vibrational dynamics of the target employed in the 
model should first be tested by comparing the calculated 
and measured angular resolved TOF spectra or energy 
transfers in either the single or multiphonon scattering 
regime, or in both if the data are available, and then 
applied to the calculations of total energy transfer. 



III. COMPARISON OF THEORETICAL RESULTS 
WITH HAS-TOF DATA 
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A. He^Cu(OOl) collisions 

The vibrational properties of th£L.GLL(001) surface have 
been extensively studied by HAaljliJ and these exper- 
iments have provided plentiful information on its vibra- 
tional dynamics whose characteristics are in many as- 
pects also relevant to other clean monocrystal metal sur- 
faces. This particularly applies to the pronounced inten- 
sities of the Rayleigh wave (RW) and the longitudinal 
resonance mode (LR) which dominate the spectra in the 
single phonon scattering regime. However, these specific 
features of the single phonon scattering spectra which can 
be clearly traced across the first surface Brillouin zone 
(SBZ) are quickly washed out as the multiphonon scat- 
tering regime is approached. Here we shall focus on the 
He^Cu(OOl) multiphonon scattering data reported in 
Ref. [ because of their amenability to the analysis of 
the energy transfer. The analyses of these measurements^ 
have shown that the multiphonon spectral intensities in 
the case of planar and atomically smooth surfaces are in 
the greatest part determined by the coupling of scattered 
He atoms to surface vibrations with polarization domi- 
nantly perpendicular to the. surface, in the present case to 
the Rayleigh phononsQotj. Hence, in our assessment of 
the energy transfer in HAS from Cu(OOl) surface we shall 
follow the description of the dynamic He-surface interac- 
tion outlined in Refs. [ with an improved modelling 
of the RW dispersion and polarization based on a dynam- 
ical matrix description used in connection with the Cu 
substratell3, and of the static He-surface potential pre- 
sented in Ref. [|l7|. 

Following the approach outlined in Sec. |l^ and Refs. [ 
|l8| we first calculate the multiphonon scattering spec- 
tra A'ki.Ts (e, AK) and compare them with the exper- 
imental HAS-TOF dataN to test the validity of the 
model potentials and vibrational dynamics parameters 
employed in the calculations. Then we use this as an in- 
put to calculate the angular resolved energy transfer (|^) 
which can again be compared with the analogous quan- 
tity obtained from expression (^) by a direct integration 
of the corresponding experimental TOF spectra and their 
first moments. 

A comparison of the theoretical and experimental an- 
gular resolved energy transfer obtained from (|^) by using 
expression (^) and the datapoints of the measured scat- 
tering spectra, respectively, is given in Fig. |l] as a func- 
tion of the projectile angle of incidence and fixed incopa- 
ing beam energy of 113 meV used in the experimentstj. 
The obtained agreement between experimental and theo- 
retical results enables the application of the thus verified 
calculated quantities (H) and in the evaluation of the 
total energy transfer using expression (^. These calcu- 
lations are reported in Sec. 



B. He^Xe/Cu(lll) collisions 

As the next typical example we consider the energy 
transfer to the adlayers of noble gas atoms adsorbed on 
metallic substrates because their vibrational properties 
have attracted considerable attention in the past two 
decades. The early interest in these systems was mainly 
due to the allegedly simple vibrational properties of such 
surfacealj which could be easier analyzed theoretically 
than the surfaces of pure metals, on the one hand, or 
of the metal compounds and ionic crystals, on the other 
hand. The later interest has been associated with the 
reduced effective dimensionality of adsorbed monolayers 
and the possibility to study phase transformations and 
thermodynamics of supported quasi-two-dimensional sys- 
tems. The most recent interest is connected with the 
development of the scanning tunneling microscope tech- 
niques which enable manipulations with adsorbates on 
the atomic scale. 

A prototype system of a rare gas adlayer is the com- 
mensurate (VS X \/3)R30^ monolayer of Xe atoms ad- 
sorbed on Cu(lll) surf ace a t substrate temperatures be- 
tween ~47 and ~80 Klijiia. The adlayer vibrations in 
this system arc almost completely decoupled from the 
vibrational dynamics of the underlying substratdlS, i.e. 
with a negligible admixture of the substrate modes and 
hence surface localized. Three Xe-adlayer induced modes 
or phonons are characterized according to their polariza- 
tion properties. The first is a shear vertical or frustrated 
translation or S-mode with polarization perpendicular to 
the surface and very flat, practically negligible disper- 
sion {ojsiQ) = 1^0 = 2.62 meV). The second is the in- 
plane dominantly longitudinally polarized L-mode which 
strongly disperses over the first SBZ of the superstructure 
and exhibits a zone center gap of the order of 0.5 meV. 
The third is the dispersive transversely polarized or shear 
horizontal (SH) mode with the in-surface-plane polariza- 
tion dominantly perpendicular to the phonon wavevec- 
tor. This mode exhibits the same zone center gap but 
its frequency lies always below that of the L-mode. The 
only deviation from the described behaviors occurs in a 
very narrow range of the phonon wavevector Q where 
the dispersion curves of these modes intersect (or rather 
exhibit an avoided crossing) with the dispersion curve of 
the RW of the Cu substrate. In this case their localiza- 
tion may extend over the first few substrate layers with 
the additional effect of mixed polarization of the modellS. 
Outside the region of avoided crossing with the substrate 
RW these modes, besides being localized in the adlayer, 
are pure in character (negligible mixed polarizations) and 
their polarization eigenvectors are mutually orthogonal. 
Hence, this system lends itself as a particularly conve- 
nient one for studying the energy transfer to pure surface 
localized vibrations. 

The results of the calculations following the general 
procedure analogous to the one employed in the preced- 
ing subsection are depicted in Fig. |[ Here the angular 
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resolved energy transfer is plotted as a function of the 
substrate temperature for several incoming He atom en- 
ergies and fixed experimental incident angle. We again 
find that the vibrations with dominant component of po- 
larization perpendicular to the surface, here the S-modes, 
give preponderant contribution to the energy transfer. 
The agreement between the thus obtained experimental 
and theoretical values is good, thereby providing a re- 
liable starting point for the calculations of total energy 
transfer reported in the next section. 

C. He^Xe(lll) collisions 

As the third example we consider the energy trans- 
fer in He atom scattering from a van der Waals crys- 
tal, viz. the (111) surface of Xe multilayer condensed 
on Pt(lll). Recent HAS angular distribution (diffrac- 
tion) and TOF studies of this systerrcj reveal an ordered 
and almost defect-free surface sustaining several types of 
well resolved vibrational modes. This enables us to apply 
the above described theoretical approach to this collision 
system as well. In particular, we were able to set up a 
dynamical matrix describing the vibrational properties of 
a Xe slab consisting of 60 layers and bounded by (111) 
surfaces by using the force constants calcpiated from the 
available Xe-Xe gas phase pair potentialg22l. The disper- 
sion curves of the various phonon modes obtained from 
such a dynamical matrix are in a fairly good agreement 
with the ones-jderived from the measured He— >Xe(lll) 
TOF spectraB. 

Application of expression j-©) to the measured 
He^Xe(lll) scattering spectraEJyields the values of the 
angular resolved energy transfer shown in Fig. |^ as a 
function of the angle of incidence and fixed incoming He 
atom beam energy of 10 meV used in the experiments. 
We then calculated the angular resolved scattering spec- 
trum (Ecu (3)) in the distorted wave Born approxima- 
tion limitQ pertinent to the scattering regime in which the 
experimental data were recorded. This was used to ob- 
tain the corresponding theoretical energy transfer from 
Eq. (2) and the result, which is presented in Fig. 3, 
is consistent with the values deduced from the experi- 
mental spectra. The general trends in the behavior of 
the energy transfer are similar to those characterizing 
the He^Cu(OOl) collision system despite the large dif- 
ferences in the incident energy of the scattered He atoms 
and the substrate temperature. 

IV. CALCULATIONS OF THE TOTAL ENERGY 
TRANSFER 

A good agreement between the experimental and the- 
oretical angular resolved energy transfer in HAS estab- 
lished in the preceding section for three prototype colli- 
sion systems enables us to employ the obtained results to 



estimate the total energy transfer in these experiments. 
Figure || illustrates the behavior of total energy trans- 
fer in the discussed systems as a function of the projec- 
tile incident energy over the interval which is relevant to 
HAS experiments as well as to the conditions of wind 
tunnel investigations of gas-surface interactions in space. 
It is noticeable that the magnitude of the energy trans- 
fer is surface specific in the studied scattering regime, 
in contrast to the cesults from standard classical accom- 
modation theoriesBu. In particular, the energy trans- 
fer is larger for surfaces with enhanced surface projected 
phonon density of low energy states with perpendicular 
polarization. Here it is largest for the Xe/Cu(lll) sur- 
face because the corresponding S-phonon exhibits per- 
pendicular polarization and Xe adlayer localization over 
almost the entire SBZ. This points to the dominant role 
which the surface modes with perpendicular polarization 
play in energy transfer in gas interactions with atomically 
smooth surfaces. 

Consistent with these properties we also find that, 
in the scattering regime spanned by the collision pa- 
rameters and energy interval of Fig. |^, the classical 
Baule fcpmula for semiquantitative estimate of the energy 
transferee applies only to the systems He^Xe/Cu(lll) 
and He— >Xe(lll). In these two cases the frequencies of 
surface vibrations of Xe atojas-icorresponding mode en- 
ergies not exceeding 4 meVll3'E2l) are much smaller than 
the inverse collision time and hence the impulsive scat- 
tering limit in which the Baule formula holds is reached. 
This is also in accord with independent testings of the 
applicability of the Baule formula to He atom scatter- 
ing from condensed multilayers of noble gas atoms on 
Si(lOO) surfaceE^. On the other hand, the surface De- 
bye fcequency of copper, w^", is much higher {Tiuj^^ ^24 
meWJ) and of the order of inverse collision time charac- 
teristic of the He^Cu(OOl) system in the studied scatter- 
ing regime. As a consequence, the impulsive scattering 
limit is not yet reached in the same energy interval which 
makes the Baule formula inapplicable. 

V. CONCLUSION 

We have presented a method for calculating the energy 
transfer in the quantum regime of scattering of atoms 
from surfaces. The general features of the results ob- 
tained for three prototype gas-surface collision systems 
recently studied by He atom scattering should be also 
relevant to a large variety of studies of free molecular 
flow past solid bodies performed in wind tunnel experi- 
ments. 
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FIG. 1. Comparison of the theoretical angular resolved en- 
ergy transfer in He^Cu(OOl) collisions (full line) calculated 
from expressions ^ and (^ with the values deduced from 
available experimental HAS-TOF spectra (open circles) using 
expression (0), given as a function of AO = 9i — Oso/'i and 
for the scattering conditions as denoted. Inset: Comparison of 
the experimental multiphonon HAS-TOF (open circles) and 
theoretical (full curve) scattering spectra for A6 correspond- 
ing to the experimental point denoted by arrow. 



FIG. 2. Comparison of the theoretical angular resolved en- 
ergy transfer in He^Xe/Cu(lll) collisions (full lines) calcu- 
lated from expressions (^ and ^ with the values deduced 
from available experimental HAS-TOF spectra (open sym- 
bols) using expression (^ , given as a function of the substrate 
temperature Ts for the scattering conditions as denoted. In- 
set: Comparison of experimental multiphonon HAS-TOF 
(open circles) and theoretical (full curve) scattering spectrum 
corresponding to the scattering conditions of the experimental 
point denoted by arrow. 



FIG. 3. Comparison of the theoretical angular resolved en- 
ergy transfer in He^Xe(lll) collisions (full line) calculated 
from expressions ^ and (^ with the values deduced from 
available experimental HAS-TOF spectra (squares) using ex- 
pression given as a function of AS — 9i — Osd/'^ and for 
the scattering conditions as denoted. Inset: Comparison of 
experimental HAS-TOF (open squares) and theoretical (full 
curve) scattering spectrum for AS corresponding to the ex- 
perimental point denoted by arrow. 



FIG. 4. Total energy transfer characteristic of 
He-»Cu(001), He-^Xe/Cu(lll) and He^Xe(lll) collisions 
calculated from expression (^ and plotted as a function of 
the incoming He beam energy for fixed incident angle and 
substrate temperature. Starred line is the energy transfer ob- 
tained from the classical Baule formula applied to He— >Xe 
binary collisions. 
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